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ABSTRACT 

We have discovered a «45 kpc Lya nebula (or Lya "blob") at z « 1.67 which exhibits strong, 
spatially-extended He ii emission and very weak C IV and C ill] emission. This is the first spatially- 
extended Lya-I-He ii emitter observed and the lowest redshift Lya blob yet found. Strong Lya 
and He iiA1640 emission in the absence of metal lines has been proposed as a unique observational 
signature of primordial galaxy formation (e.g., from gravitational cooling radiation or Population III 
star formation) , but no convincing examples of spatially-extended Lya-I-He ii emitters have surfaced 
either in Lya-emitting galaxy surveys at high redshifts {z > 4) or in studies of Lya nebulae at lower 
redshifts. From comparisons with photoionization models, we find that the observed line ratios in this 
nebula are consistent with low metallicity gas {Z < 10^^ — 10~^Zq), but that this conclusion depends 
on the unknown ionization parameter of the system. The large He ii equivalent width (w37±10A) and 
the large He ii/Lya ratio (0.12±0.04) suggest that the cloud is being illuminated by a hard ionizing 
continuum, either an AGN or very low metallicity stars, or perhaps powered by gravitational cooling 
radiation. Thus far there is no obvious sign of a powerful AGN in or near the system, so in order 
to power the nebula while remaining hidden from view even in the mid-infrared, the AGN would 
need to be heavily obscured. Despite the strong Lya-f He ii emission, it is not yet clear what is the 
dominant power source for this nebula. The system therefore serves as an instructive example of how 
the complexities of true astrophysical sources will complicate matters when attempting to use a strong 
Lya-|-He II signature as a unique tracer of primordial galaxy formation. 
Subject headings: galaxies: formation — galaxies: high-redshift — galaxies: surveys 



1. INTRODUCTION 

Understanding primordial galaxy formation is a ma- 
jor science driver for the next generation of large space- 
and ground-based telescopes and has inspired a substan- 
tial amount of theoretical literature due to the potential 
contribution of the first generations of stars to reion- 
ization and the early stages of galaxy evolution. Ac- 
cording to theoretical predictions, the observational sign- 
post of primordial galaxy formation is the presence of 
strong LyaA1216 and He iiA1640 emission lines, either 
due to photoionization by very low metallicity and Pop- 
ulation HI (Pop HI, i.e., z ero metallicity) star forma- 
tion (e.g.. iTumlinson et all ,2001: Schacrcr 2003, 2008) 
or due to gas cooling during gravitatio nal collapse (e.g., 
iHaiman et all 120001 : lYang et al.l l2006f) . In the case of 
Pop III stars, the strong Lya and He ii is a direct conse- 
quence of the low metallicities, where, in the absence of 
metals, H and He become the dominant line coolants for 
the gas, and of the hot effective temperatures of Pop HI 
stellar clusters, which are predicted to show a top-heavy 
IMF and low stell ar atmospheric opacity d ue to the lack 
ofmetals (e.g., lEzer fc CameronI 119711 : iBromm et"al] 
[200l . In the case of gravitational cooling radiation, Lya 
and He ii are the primary ways for pristine gas to cool 
as it is collisionally excited during gravitational collapse; 
the p redicted He ii/Lyg ratios may be as high as 10% 
(e.g., IHaiman et aLll2000l : lYalig et al.ll2006l ). Strong Lya 
and He il emission lines are commonly seen in other as- 
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trophysical sources with hard ionizing continua, but at 
normal metallicities they are generally accompanied by 
strong metal lines such as C ill] and C iv, as seen for ex- 
ample, in AGN, radio galaxy halos, Wolf-Rayet galaxies, 
or cases of sho ck ioniz ation (e.g ., Rculand et al. 200j; 
iLeitherer etliLl Il996l: iDopita fc Sutherland! Il996l) . The 
presence of strong Lya and He ii emission in the absence 
of strong metal lines has been put forward as a poten- 
tially unique observational signature of primordial galaxy 
formation. 

Although theoretical studies suggest that Pop III and 
very low metallicity star formation may persist down 
to lower redshifts, this depends on the feedback effi- 
ciency, i.e., the ability of a Pop HI stellar population 
to pollute the large-sc ale surroundings with metals (e.g., 
iTornatore et al.l|2007f ). Searches for Pop HI stars have 
understandably pushed to higher redshift (z > 4) , where 
the Pop HI star formation rate density should increase 
dram atically relative to that found in the local Uni- 
verse (IScannapieco et aLll2003[ ). Thus far, no unambigu- 
ous case of a Pop HI stellar population has been ob- 
served. Several Lya-emitting galaxy studies at z ^ 4 — 5 
have uncovered sources with unusually high Lya equiv- 
alent widths {Wrest > 240 A) - larger than that ex- 
pected from a normal stellar population - suggesting ei- 
ther a top-heav y IMF, a very low rn etallicity , and/or a 
very young age jMalhotra fc Rhoadsir2002; Rhoads et all 
l2003l ). However, the complicated radiative transfer of 
Lya in a clur npy ISM could also be r esponsible for boost- 
ing the Lya (|Finkelstein et al1l2008f ). and in these stud- 
ies, no corroborating evidence for the Pop III scenario in 
the form of a strong He ii detection was found in either 
the individual or stacked LAE spectra, leaving the mat- 
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ter unresolved (|Dawsoii et all |2004 lOuchi et al.ll2008D . 
Deep spectroscopic observations of a strong Lya-emitting 
galax y at z ^ 6.33 sho wed no evidence for He ii emis- 
sion ( Nagao et al.l 120051 ). and a more recent Lya+He ii 
dual emitte r survey at z ^ 4 — 5 found no convincing 
candidates (jNagao et al.ll2008[ ). The hmits from each of 
these studies suggest that Pop III star formation does 
not dominate in these z ^ 4 — 5 samples, prompting the 
authors to encourage searches at ever higher redshifts 
(z > 7). Additio nal high red shift Ly a+He ii su rveys are 
underway (e.g., Idi SeregoAl ighicri et al.ll2008D . 

In this paper we report on the discovery of a z « 1.67 
Lya nebula with strong, spatially-extended He ii emis- 
sion and very weak C IV and C ill] emission. This is 
the first spatially-extended source that resembles the pre- 
dicted Lya-|-He ii signature of primordial galaxy forma- 
tion. However, the system is more complex than it first 
appears. The observed line ratios suggest that the neb- 
ula may contain low metallicity gas, but this depends on 
the unknown ionization parameter of the system. De- 
tailed analysis of the spectra along with extensive multi- 
wavelength data reveals that the source of ionization is 
uncertain: the nebula is either an H ii region ionized 
by a hard spectrum source, i.e., an AGN or a very low 
metallicity stellar population, or a gravitationally cool- 
ing cloud. The fact that multiwavelength follow-up ob- 
servations are required in order to better constrain the 
source(s) of ionization and metallicity of the nebula has 
implications for Lya-|-He II searches at higher redshift. 
In Section 2, we summarize the systematic search for Lya 
nebulae that led to this discovery and our observations 
and reductions. Section 3 contains a discussion of the 
observational results, and Section 4 details our analysis 
of the physical properties of the nebula. In Section 5 we 
discuss the implications of this discovery for the ongoing 
high redshift Lya+He ii surveys, and we summarize our 
conclusions in Section 6. 

We assume the standard ACDM cosmology (f2M=0.3, 
riA=0.7, h=0.7); 1" corresponds to a physical scale 
of 8.47 kpc at z = 1.671. All magnitudes are in 
the AB system. Unless otherwise stated. He ii refers 
to He iiA1640A, C iv to C ivAA1549,1550, C iii] to 
C iii]A1909, and Ne iv] to Ne iv]A2424. 

2. OBSERVATIONS & REDUCTIONS 

In this section we discuss the observations that led to 
the discovery of this Lya-f He ii nebula and the methods 
used to process the imaging and spectroscopic data. 

2.1. The Search 

Large Lya nebulae (or Lya "blobs") — large 
(~100 kpc) clouds of gas emitting strongly in Lya — 
are thought to be sites of ongoing galaxy formation 
and have been found in small numbers, primarily at 
2 < z < 3. Early theoretical work suggested that 
these nebu lae could be exampl es of gravitationally cool- 
ing clouds ijHaiman et al.l l2000^. and more recent cosmo- 
logical hydrodynamic simulations indicated that cooling 
clouds should be detectable as Lya-I-He li nebulae, al- 
though the specific predictions proved uncertain due to 
the treatment of star-forming gas (Yang et al. 2006; 
Y. Yang 2008, private communication). Several of the 
largest Lya ne bulae, including o ne with strong He ii and 
C IV emission (|Dev et al.ll2005l) . have since been shown 



to be powered instead by AGN, spatial ly-extended star 
formation, or some combination (e.g.. |Dev et al.l l2005t 



iMatsuda et"al] 120071: iGeach et al] 120071 1 Two groups 
claim to have discovered Lya nebul ae that are powered 
by gravitational coohng radiatio n (Nilss on et al.l l2006t 
ISmith fc Jarvisl2007HsInith et al.>2008. ). but neither case 
shows strong He ii emission. 

Lya nebulae are extremely rare objects and have of- 
ten been found using deep narrow-band imaging sur- 
veys of known galaxy overdensities. As such, their space 
density, particularly at the bright end of the luminos- 
ity function, is largely unconstrained. In order to carry 
out an efhcient but unbiased survey of a large cosmic vol- 
ume, we have designed a systematic morphological search 
for spatially-extended Lya nebulae using the broad-band 
imagin g from the NOAO Dee p Wide-Field Survey (ND- 
WFS; I Jannuzi fc Devi fl999l) Bootes Field. This field 
has been imaged in Bw, and /-band with median 
5a point-source depths of w27.1, 26.3, and 25.8 mag 
(AB), respectively. In the sub- field relevant to this pa- 
per, the Bw, R, and /-band 5a point-source depths are 
27.5, 26.1, and 26.0 mag for 2.3, 1.7, and 2.8 hrs of inte- 
gration, respectively. Candidates were selected from the 
broad-band data using a morphological and color selec- 
tion algorithm, and spectroscopic follow-up was used to 
rule out low-redshift interlopers. The advantage of us- 
ing broad-band data to search for line-emitting sources, a 
seemingly crude approach, is the enormous comoving vol- 
umes (-10^ hfo^ Mpc3) that can be surveyed efficiently 
using publicly available data over wide fields. The suc- 
cess of the present survey relied heavily on the depth of 
the broad-band NDWFS imaging and the darkness of the 
sky within the /?vi^-band, against which strong line emis- 
sion can dominate the flux even within the very broad 
Bw filter. 

A full discussion of the search algorithm and results 
will be discussed elsewhere (Prescott et al. 2009, in 
preparation). Here we report on the discovery of a new 
Lya-fHe ii nebula at z ss 1.67 (hereafter denoted PRGl) 
located in the Bootes field at 14:35:12.439 -^35:11:07.16 
(J2000). Unlike the other Lya sources found in our sys- 
tematic search, the Lya-I-He ii nebula presented here 
was selected by the search algorithm not because of 
its Lya emission, which lies outside the Bw filter, but 
instead because of the strong and diffuse blue contin- 
uum emission («92%) and spatially-extended He ii emis- 
sion («8%) within the bandpass . Postage stamps from 
GALEX (FUV & NUV: iMarti^et al...200& ). the NOAO 
Deep Wide-Field Survey {Bw, R, and /), and the Spitzer 
Deep Wide-Field Survey (SDWFS; IRAC 4.5^m, 
5.8/im, and 8.0/im; Ashby et al. 2009, submitted) are 
shown in Figure [T] The GALEX limits are 0.36 /iJy in 
the NUV and FUV bands; the SDWFS limits are 3.2, 
4.4, 25.5, and 25 /iJy {5a). The MIPS coverage of this 
region shows no detection with a la rms limit of 51 /iJy 
(E. Le Floc'h 2008, private communication). 

2.2. Spectroscopic Follow-up 

We obtained spectroscopic follow-up observations us- 
ing the MMT and the Blue Channel Spectrograph dur- 
ing UT 2008 June 8-9 as part of our systematic search for 
Lya nebulae. We used a 1.5x120" (unvignetted) slit and 
the 300 1/mm grating (Ac « 5713A, AA « 3100-8320A). 
We chose a slit orientation that spanned the longest 



A Lya+Heii Nebula at z w 1.67 



3 



dimension of the diffuse emission (PA=81.2°, observed 
near transit), as shown in Figure [U while also intersect- 
ing a nearby bright object. We moved the target along 
the slit by ~5" inbetween exposures. 

Conditions during the first night were clear and stable 
with I'.'O seeing; the second night was clear but with vari- 
able seeing (>1'.'3) and high winds which caused shaking 
of the telescope pointing. We show in Section 3.2 that 
due to a slight pointing offset and wind-shake, the data 
from the second night sample a different spatial region 
within the nebula and are contaminated at some level by 
sources nominally off the slit, in particular a red com- 
pact source to the WNW (Source A in Figure [1]). For 
this reason, data from the two nights were reduced and 
analyzed separately. The most robust Lya and He ii flux 
measurements are from the first night, when conditions 
were excellent during the single half hour exposure. An 
additional 1.5 hours of integration were obtained during 
the second night. The Lya flux measurement from Night 
2 shows a 25% loss relative to that from Night 1. In the 
remainder of the paper we use the Lya and He ii flux 
measurements from Night 1 but include the line ratios 
derived from both nights. 

The data were reduced using IRAF'^ . After performing 
the overscan and bias subtraction, we corrected flatfield 
exposures for the response of the internal 'Bright Con- 
tinuum' flatfield lamp by dividing out the median along 
columns and then applied the flatfield correction. We 
used twilight flats to determine the illumination correc- 
tion for the science frames. Cosmic rays were removed 
using xzap^ . The wavelength solution was determined us- 
ing HeArNe and HgCd comparison lamps, with an rms 
of «0.17A. We corrected the data for a slight systematic 
offset in the night sky lines; the night sky line wave- 
lengths in the final spectra are correct to ±0.3A. Flux 
calibration was based on observations of the standard 
stars BD-l-33 2642 and Wolf 1346^ We applied a grey 
shift (<0.08 mag) and fit the sensitivity function using 
extra care at the blue end of the spectrum because the 
Lya line at «3250A hes only 56.3 pixels (109. TA) from 
the edge of the chip. The instrumental resolution mea- 
sured from the Hg iA4047 line is 3.6A, and the tilt within 
the aperture is < ±5.24 km s~^ over the region of the 
nebula. 

3. RESULTS 

3.1. Lya and He II Emission 

The final 2D and ID spectra from the first and second 
night are shown in Figures [2][3l The spectra show strong 
Lya and He ii emission lines, both of which exhibit sim- 
ilar kinematic structure in the 2D spectrum. The mea- 
sured fluxes and flux limits were derived separately from 
each night's data using a 1.5x5.0" aperture, chosen to 
maximize the signal-to-noise ratio of the He ii measure- 
ment (Table [T]) . Faint continuum emission is detected 
in the spectra (Figure [3|) . We scaled the spectra by a 
factor of 1.4 and 1.6, respectively, in order to match the 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

* http;/ /iraf.noao.edu/iraf/ftp/iraf/extern/xdimsum020627 

5 KPNO IRS Standard Star Manual 



continuum fluxes measured from the much deeper ND- 
WFS broad-band data within the region covered by the 
slit (Tabled]). 

The Lya and He ii luminosities from Night 1 are 

Lbya = Lhya.ap X fgeo X /profile ~ 5.4 X lO'^^ Crg S"^ and 
LhcII = LhcUmp X fgeo « 4.0 X 10"*^ erg S~^ LLya.ap = 

9.3 X lO'^^ erg s^^ and LneiLap — 1-2 x lO''^ erg s~^ are 
the luminosities measured within the spectroscopic aper- 
ture, fgeo is the geometric correction factor between the 
spectroscopic aperture and full extent of the nebula, and 
f prof He corrects for asymmetry in the Lya profile due to 
blue side absorption. We discuss the estimation of these 
correction factors below. 

The geometric correction factor fg^o was derived from 
a comparison of the spatial extent of the nebula measured 
in the Lya line versus the extent in the Bw imaging. In 
our best seeing Night 1 data, where we have accurate 
spatial information and low slit losses, the Lya is ex- 
tended by «5" (42.3 kpc) with a fairly sharp truncation 
at large radii, but with a possible extension towards the 
west (right of center. Figure S]). The He ii emission is 
weaker than the Lya by a factor of 10, but it appears 
that at lower signal-to-noise the He ii nebula extends to 
roughly the same radius as the core of the Lya nebula 
(no westward extension). The Night 2 data show a Lya 
extent of «6" and a He ii extent of «5" but suffer from 
slit losses and degraded spatial resolution due to poor 
seeing and wind-shake. The extent of the Lya nebula 
in the spectroscopic observations is in rough agreement 
with the spatial extent of the diffuse emission measured 
from the broad-band Bw data along the position of the 
slit («6.6", «56 kpc above a Bw surface brightness of 
4.5x10"^^ erg s~^ cm~^ arcsec"^). The full area of the 
diffuse emission in the Bw imaging is « 26 square arcsec- 
onds. Assuming that the Lya emission is distributed sim- 
ilarly to the Bw continuum emission, we estimated that 
a geometric correction factor of fgeo = 3.4 is required to 
obtain the total Lya flux from the nebula. This correc- 
tion is very approximate; narrow-band imaging and/or 
spatially-resolved spectroscopy will be required to accu- 
rately account for the contributions of line and contin- 
uum emission to different portions of the nebula. 

We derived the Lya blue side absorption correction 
fprofiie using a comparison of the Lya and He ii line 
centroids. The Lya profile is fairly symmetric (Figure [5]), 
but the center of the Lya line is offset to the red from 
the systemic redshift, as determined from the centroid 
of the He ii emission line, likely due to absorption. We 
estimated the amount of blue side absorption of Lya by 
mirroring the red side of the Lya profile across the line 
centroid, and derived a factor of 1.7 (Night 1) and 1.8 
(Night 2) increase in the Lya flux. We therefore used 
a correction factor of f profile = 1.7 to obtain the final 
Night 1 Lya luminosity. 

The nebula shows clear velocity structure in both the 
Lya and He ii lines. The lines are resolved, showing Lya 
and He ii velocity dispersions of cr^ ~ 400 and 250 km 
s~^, respectively, corrected for the instrumental resolu- 
tion. Figure [5] shows the velocity profile of the Lya line 
from Night 1 and Night 2 derived using 2 pixel (0.56") 
extractions. The profile is relatively smooth but flattens 
on the east (left) of center. The spatially-resolved Lya 
velocity dispersion is essentially constant across the neb- 
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ula. 

3.2. C IV, C III], and Ne IV] Emission 

Despite the strong Lya and He ii emission seen in 
the Night 1 data, there is no detection of C ill], C IV, 
or Ne iv]. The la upper hmits on the hne ratios are 
C iv/He IK 0.23, C iii]/He ii< 0.19, and Ne iv]/He ii< 
0.23. Due to the excellent and stable observing con- 
ditions during Night 1, the Night 1 spectrum provides 
the most accurate flux measurements for Lya and He ii 
along with higher resolution kinematic information and 
the strongest limits on the C ill], C iv, and Ne iv] emis- 
sion from the source. In contrast, the Night 2 spec- 
trum shows weak C iii] and marginal C iv and Ne iv] 
emission at the same redshift. The resulting line ra- 
tios are C iv/He ii = 0.36, C iii]/He ii < 0.82, and 
Ne iv]/He ii < 0.49. The Night 2 data provide addi- 
tional constraints on the line ratios but must be treated 
with care due to the Night 2 observing conditions. Due 
to the poor seeing and the wind-shake of the telescope, 
the Night 2 spectrum suffers from slit losses as well 
as contamination from nearby sources nominally off the 
slit, most importantly from a compact red source to the 
WNW of the target center (Source A; see Figure [1]). 

From a comparison of the line ratios derived from each 
night, we argue that the metal line emission is not from 
the same spatial location as the Lya and He ii, and may 
instead be associated with the region closer to Source 
A. On Night 2, the C iii]/He ii ratio was 0.82; thus if 
the emission were perfectly cospatial, we should have de- 
tected C III] on Night 1 at «5.1xl0~^'' erg s~^ cm~^, 
roughly the same significance as He li. Instead, we can 
rule out C III] emission at the 4.3cr level. For Ne iv], 
we should have detected it at 3.0x10^^^ erg s^^ cm^^ 
on Night 1, which disagrees at the 2.2a level with our 
Night 1 result. We cannot make a robust comparison for 
C IV, as it is only detected at the 2a level even on Night 
2, but we make the assumption that all the metal line 
emission originates from the same source. In contrast, 
the He ii/Lya line ratios are consistent between the two 
nights (0.12±0.04 and 0.13±0.02, respectively), indicat- 
ing that the Lya and He ii are indeed cospatial within 
the region sampled by these observations even though 
the metal line emission varies spatially. 

The offsets in the spatial profiles of the lines are con- 
sistent with the idea that the data from the two nights 
sampled different spatial regions. If the C iv, C ill], and 
Ne iv] lines are from the region around Source A and if 
the wind-shake was consistently perpendicular to the slit, 
there should be a «5.2 pixel (1.5" West) offset between 
the spatial centroid of these lines and that of Lya. How- 
ever, this offset will vary by an estimated ~ ±3.6 pixels 
(Ril") or more depending on the direction of the tele- 
scope wind-shake relative to the angle of the slit during 
a given exposure. We do see a difference between the 
two nights when we look at the spatial profiles extracted 
in 2 pixel spatial bins along the spatial direction (Fig- 
ure [4]) ; the Night 1 profile is skewed towards the East 
(left, away from Source A), whereas the Night 2 profile 
is peaked closer to the position of Source A, suggesting 
contamination. We use these spatial profiles to compute 
flux-weighted mean spatial centroids for each emission 
line separately. The He ii position is consistent with that 
of Lya to within 1 a (Ax=1.41±1.54 pixels, 0.40±0.43"). 



The C III] offset is 4.03±2.07 pixels (1.13±0.58") in the 
direction of Source A. (Due to the intrinsic faintness of 
the C IV and Ne iv] lines, the computed offsets are not 
statistically significant.) Follow-up observations will be 
required to resolve this issue, but it appears that the 
C III] line is offset from the spatial centroid of the Lya 
at the 1.9ct level, in the right direction and at roughly the 
correct position to be explained by contamination from 
the region closer to Source A. 

Given the pointing uncertainty and inevitablity of con- 
tamination from sources off the slit due to the unstable 
conditions on Night 2, the discrepancy in the spatial pro- 
files between the two nights, the evidence that the C ill] 
line shows a spatial offset consistent with the position 
of Source A, and the fact that the Night 1 data are in- 
consistent with the fluxes of Ne iv] and C ill] measured 
on Night 2 at ^2-4f7, we argue that much of the metal 
line emission derives from the region around Source A. 
We cannot rule out that some is emitted further out in 
the nebula, but even if it does, it arises from a region 
spatially distinct from the region observed on Night 1. 
Combining the data from the two nights would not be 
appropriate due to the different spatial sampling of the 
observations and these intrinsic spatial inhomogeneities. 
In the absence of more definitive data, we use the metal 
emission line measurements from both nights as inde- 
pendent upper limits on the emission coming from two 
possibly distinct regions of the nebula. 

4. DISCUSSION 

In this section we derive estimates for the physical 
properties of PRGl and use photoionization models to 
gain insight into the gas metallicity and the possible 
source(s) of ionization. Despite the strong Lya-I-He li 
signature and weak metal line emission often associated 
with primordial phenomena, we find that it is not possi- 
ble to make an unambiguous determination of the source 
of ionization. The line ratios are consistent with a nebula 
comprised of low (but non-zero) metallicity gas, irrespec- 
tive of the nature of the ionizing source, but this conclu- 
sion depends on the unknown ionization parameter of the 
system. 

4.1. Physical Properties of the Nebula 

Our discovery data can be used to put constraints on 
the physical properties of the nebula. Using the emission 
measure of Lya, we can estimate the electron density: 

^^v"' = if^-^P^'^f'^h^Hpail (1) 
« l.2^nlfVhvHpa^il (2) 

where jhya and jup are the emission coefficients for 
Lya and H/3, respectively, rip and are the proton 
and electron number densities with rtg « 1.2np (the 
factor of 1.2 accounts for the contribution of electrons 
from doubly-ionized Helium), / is the volume filling fac- 
tor of the nebula, V is the volume of the nebula, h 
is Planck's constant, vnp is the frequency of if/3, and 

effective recombination coefficient for Hj3 
(jOsterbrockl I1989D . We approximate the nebula as a 
sphere with radius R ~ 28.0 kpc. The Lya luminos- 
ity extrapolated to the entire nebula and corrected for 
blue side absorption (5.4 x 10^^ erg s~^) corresponds to 
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an electron number density of rie ~ 0.094/" 
29.7(//10^^)^°'^ cm^'^, where we have used a typi- 
cal value for / derived for the line-emi tting regions in 
cluster cooling flows 10~^; e.g., iHeckman et al.l 
Il989t ). This corresponds to an ionized gas mass of 
M,on=l.25mpnJV =8.4x 108(//10~5)0.5 Mq. 

Similarly, the He ii emission measure can be used to 
estimate the He"*""*" and electron densities, assuming the 
cosmic mass fraction of He: 



^ J A 1640 



(14)n^ 



A4686 
eff 



A4686 



(3) 

(4) 



where jAi640 and Jmssg are the emission coefficients, 
njjf,++ and Ue are the He"*""*" and electron number den- 
sities with Tie « 1.2np « 14n//g++, ;^a4686 is the fre- 
quency of He IIA4686, and at^ Lfi is the He iiA46 86 ef- 
fective recombination coefficient (|Osterbrockl ll989^ . The 
He II luminosity extrapolated to the whole nebula (4.0 x 
10*^ erg s^^) corresponds to a He++ number density of 
'^He++ ~ 1.0(//10^^)~°'^ cm""^ and an electron num- 
ber density of ~ 14.5(//10^''')^'^'^ cm^'^. This corre- 
sponds to an ionized gas mass of A/io„=1.25mpne/T^ = 
4.1x 10®(//10~^)"-^ M0, roughly consistent with the Lya 
estimate. 

If we make the naive assumption that the veloc- 
ity spread results from rotation with Vc « = 
350 km s~^ at a radius of 28 kpc (3.3"), we 
estimate the mass interior to this radius to be 
^Irot =8.0x lO-'^^sin^z Af©, where i is the inclination of 
the system. If the velocity dispersion is due to random 
motions of small clouds within the system, we estimate 
the dynamical mass (using the more kinematically robust 
He II line; = 245 km s~^) to be Mrand = 1-9 x 10^^ 
Mq. 

The spatial extent of the Lya (H+) and He ii (He++) 
emitting regions appear to be similar in size, certainly 
within a factor of two {Rh+ / RHe++ S 2)- This is some- 
what surprising when we consider a simple Stromgren 
sphere argument. Assuming a hard central ionizing 
source (an AGN or a Pop III SED) embedded within an 
infinite spherical H-|-He cloud, the predicted size ratio of 
the H"*" and He+^-emitting regions is Rh+ / RHe++ ~ 9 — 
17, an order of magnitude higher than we observe. This 
suggests that either the H+ region is density-bounded, 
that the sources of ionization are distributed throughout 
the nebula, or that the assumption of spherical symmetry 
is invalid (e.g., the cloud is illuminated from the outside). 

The observed Lya luminosity {LLya — 5.4x10'*'^ 
erg s""'^) corresponds to a H- ionizing photon flux 
(13.6 eV < E.^ < 54.4 eV) of: 



Q{H) 



L 



Lya 



1 



4.9 X 10'''' photons s~ 



hvLya 0.68 



(5) 



where we have assumed that the fra ction of H io nizing 
photons conver ted into Lya is 0.68 (|Spitzeij[l978l) . We 
note that this is likely a lower limit due to a number 
of considerations: Lya is highly susceptible to resonant 
scattering and is easily destroyed by dust, and Lya for 
this system is observed at «3250A, a wavelength regime 
that suffers from very low atmospheric transimission and 
poor CCD sensitivity, making accurate flux calibration 
difficult. 



From the observed He ii emission [LhcII = 4.0x10''^ 
erg s~^), we calculate a He+-ionizing photon flux (E-y > 
54.4 eV) of: 

Q{He+) = ~ 6.2 X 10^3 p/jotons s-\6) 



•"Hell 



where a^/jj — 



1.53 X IQ-^^ cm3 s^i (case B; 



100 cm-3, 10* K: lStorev fc Hummei]ll995h and a^¥° 



,,4686 jAi640 i^xiese 



8.08 X 10-13 s-i (case B; 



IOsterbrocklll989l) . 

The large value of Q(He"'")/Q(H)=0.13 is strong evi- 
dence that the source is illuminated by a hard ionizing 
continuum. The prediction for a Pop II stellar popu- 
lation (insta ntaneous burst , Salpeter IMF, 1-100 Mq, 
Z = 0.001: ISchaereJ [20031) is Q(He+)/Q(H)= 0.0004; 
while our Lya measurement may be compromised by 
radiative transfer effects and poor flux calibration at 
the edge of the atmospheric cut-off, it would need to 
have been underestimated by more than a factor of 
100 in order for the observed Q(He"*')/Q(H) to match 
that of a normal Pop II stellar population. The ob- 
served Q(He+)/Q(H) is in the range expected for AGN 
or Pop HI stars. A t ypical AGN temp late has a ratio of 
Q(He+)/Q(H)-0.09 (|Elvis et al.l 11994 V more consistent 
with the observational constraint. While normal metal- 
licity stellar populations will have very little flux above 
the He"*" ionization edge. Pop HI stars are predicted to 
have much harder spectra due to high temperatures, low 
stellar atmospheric opacity, and a top-heavy IMF. Strong 
mass loss from these stars would cause higher effective 
tempertu res, boos t ing th e hard ionizing photons even 
further. ISchaereil (|2002f ) calculated a suite of Pop III 
models, both with and without strong mass loss. For 
individual high mass stars (M>80-300 Mq) they predict 
Q(He+)/Q(H)>0.022-0.11. Models with mass loss yields 
ratios of Q(He+)/Q(H)>0.06-0.17 (80-300 Mq). How- 
ever, such large Q(He+)/Q(H) ratios persists for only a 
few Myrs for instantaneous burst models. Constant star- 
forming Pop III models (with no mass loss) integrated 
over a range of IMFs yield Q(He+)/Q(H)< 04, and 
the v alue decreases with increasing metallicity (jSchaered 
[2003h . There are numerous uncertainties in these esti- 
mates, but broadly speaking the observed Q(He"'")/Q(H) 
ratio is in the range populated only by AGN and the very 
lowest metallicity stellar populations. 

4.2. Photoionization Modeling 

The ubiquity of AGN, the association of AGN and ex- 
tended emission line regions (EELRs), and the presence 
of spatially varying metal line emission suggest that an 
AGN is a plausible source of ionization. At the same 
time, the observed line ratios are highly unusual and 
have been tied in the theoretical literature to primor- 
dial galaxy formation processes (the presence of Pop III 
stars or gravitational cooling radiation). A compari- 
son of the line ratios from PRGl with those of radio 
galaxy EELRs and other Lya nebulae is shown in Ta- 
ble [3l While the Lya/He ii ratios are comparable, other 
than the C ill] /He ii ratio on Night 2, all the metal line 
ratios from this source are at the low end or lower than 
the range seen in EELRs. Furthermore, this source shows 
very different line ratios than those seen in another large 
radio-quiet Lya nebula at z w 2.7 (jPev et al.ll2005f ). 
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In order to gain greater insight into the metaUicity 
and possible ionization sources for the nebula, we used 
Cloudy^ to model simple, constant density gas clouds be- 
ing illuminated by AGN, Pop III, and Pop II SEDs and 
predict the resulting line ratios and conti nuum emission. 
The A GN template is taken from Mathews &: Ferlandl 
(|T98l {F^ cx at 26 < Iw < 56 eV, oc v-^ ° 

at 56 < ht^ < 3 65 eV ), and the Pop III spectra are 
iTumlinson eTall (|2006f ) 1 Myr models {Z = 0, a top- 
heavy IMF peaked around lOM© with (t=1.0, i.e., their 
case A). The Pop II case is a 1 Myr, Z = 0.001, 
Sal peter IMF, instantane ous burst model from Starburst 
99 (|Leitherer erallll996l ). 

The strong, spatially-extended Lya-I-He ii emission 
and weak, spatially-variable C ill] and C iv emission in 
our discovery spectra and the blue, spatially-extended 
continuum emission measured from deep broad-band 
imaging provide constraints on the metaUicity of the gas 
and on the slope of the ionizing continuum in the system. 
He+ and C+'^ have similar ionization potentials (54.4 eV 
and 47.9 eV, respectively), so a comparison of He ii and 
C IV puts constraints on the metaUicity that are less de- 
pendent on the ionizing continuum slope. The Lya and 
C III] emission lines (relevant ionization potentials of H 
and C'^'^: 13.6 eV and 24.4 eV) provide additional con- 
straints on the slope of the ionizing continuum. 

The observed lines will also depend on the ionization 
parameter of the system. A very rough estimate based on 
our discovery data is U = (j){H)/{nHc) > 0.0002, where 
(/)(H)=Q(H) / {ttR'^) is the surface flux of ionizing photons, 
hh is the total hydrogen gas density, and R is the radius 
of the cloud. We have assumed our observed quantities: 
Q(H)=4.9xl05'' photons s-\ njj - rie - 29.7 cm^^, 
and i? « 28 kpc. However, this estimate is uncertain 
and most likely a lower limit. We have argued that the 
Lya flux is likely underestimated and raised the possi- 
bility that the system is density-bounded, both of which 
will cause an underestimate of the ionizing photon flux 
(possibly by an order of magnitude). In addition, the 
geometry of the system is clearly complex, and the esti- 
mated density depends critically on the assumed value of 
the filling factor. The density estimates in Section 4.1 are 
larger than typical values for the ISM {uh ~1 cm~^); if 
the cloud is in fact more similar to typical ISM densities, 
the ionization parameter would increase by over an or- 
der of magnitude. A more sophisticated treatment of the 
ionization parameter is beyond the scope of this analysis, 
so for the purposes of this paper we modeled a reasonable 
range of ionization parameters: Log U=[-3, -1, 0]. 

4.3. MetaUicity of the Gas 

In order to explore the range of parameter space al- 
lowed by our observations, in Figure [7] we plot line diag- 
nostics for models with a range of ionization parameters, 
metallicities, and ionizing SEDs, comparing them to our 
observed limits on the line ratios of the nebulosity from 
Night 1 and Night 2. 

For low ionization parameters (Log U < — 1), it is pos- 
sible to produce the observed line ratios with an AGN 
SED illuminating solar metaUicity gas (Figure [71 blue 
plus signs). For higher ionization parameter (Log U > 

^ Calculations were performed with version 07.02.02 of Cloudy, 
last described bv lFerland et al.l ||T998|) . 



— 1), the observed ratios require low metaUicity gas. In 
the case of an AGN SED, the C iv/He ii and C iii]/He ii 
Umits imply Z < IO^^Zq. Models with Pop HI ioniza- 
tion lead to lower metaUicity estimates of ^ < lO"'^^©. 
(For the Pop II model even lower metallicities would 
be required, but this case is highly unlikely given the 
large Q(He+)/Q(H) ratio, as discussed in Section 4.1.) 
To put this into context, these metaUicity estimates are 
at or below the lowest limits for we ak Mg ii absorbers 
at ^4 < z < 2.4 (>10-2-10-2-5; iLvnch fc Charitoiil 
|2007( ). Studies of the most metal-poor damped Lya ab- 
sorbers have shown t hat none have met allicities lower 
than [Fe/H] > -3 (jPettini et al.ll2008D . Lya forest 
clouds at z w 3 with N(Hi)>10^^ cm~^ are uniformly 
metal-enriched with carbon abundances of ~ 10~^ Zq, 
and observations of lower column density Lya forest 
clouds indicate that there may be universal minimum 
metaUicity of 3.5 xlQ-"^ Z© th at is roughly constant 
from z « 2 - 6 (jSongailal [200l . If the metaUicity of 
PRGl is similar to the lowest metallicities measured from 
absorption-line studies, this would be the first time such 
a system has been seen in Lya and He ii emission. 

While the weak metal line emission suggests a low 
metaUicity system, the metaUicity estimates are uncer- 
tain due to a variety of factors, none of which are 
well-constrained by the current data, e.g., the ioniza- 
tion parameter, the geometry of the cloud and ionizing 
source(s), and the intensity and spectrum of the source(s) 
of ionization. Putting stronger constraints on the metal- 
licity of the nebula will require deeper spectroscopy and 
more detailed photoionization modeling, which is beyond 
the scope of the current paper. 

4.4. Source of Ionization 

The large He ii equivalent width and large He ii/Lya 
ratio is strong evidence that the nebula is powered by a 
hard ionizing continuum. We discuss each of the possible 
ionization sources in turn: AGN, shock ionization, Wolf- 
Rayet stars, low metaUicity (Pop II) and zero metaUicity 
(Pop HI) star formation, and gravitational cooling radi- 
ation. 

4.4.1. AGN 

An AGN can produce high He ii/Lya ratios and weaker 
C IV and C ill] emission lines, particularly if it is il- 
luminating a low metaUicity cloud. This scenario is 
certainly plausible, as emission line halos (e.g., Lya, 
C III], He II, C iv) around radio galaxies have been 
known for some time, arising predominantly from a 
combination of jet-ISM interactions and scattered light 
from the AGN or from spatially-extended star forma- 
tion (e.g. McCar thv et all Il98'7l: Iva n Oiik e t all 119961: 
Dev et al.iri997. : .Villar-Martin et al.ll2003: Reuland et all 
20031) . Unlike our Lya-I-He ii nebula, however, these 
gaseous haloes are predominantly metal -enriched, with 
strong C III] arid C ly emission (e.g., Reul and et al.ll2007l : 
Maxfield ctiDIlOOl). 

While there is no compact source visible in the center 
of the nebula, there are several compact sources in the 
vicinity of the nebula that could in principle harbor an 
AGN (Figure [T]). For two reasons we believe that if there 
is an AGN in the system, it must be at the position of 
Source A. First, we find that even if all the nearby sources 
were AGN, Source A would contribute the vast majority 
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of the ionizing photon flux. When we sca le the standard 
quasar template from lElvis et"al] (|1994l ) to match the 
measured Bw flux from each source and calculate the 
corresponding ionizing photon flux, the net contribution 
is only 2% of the ionizing photon flux required to explain 
the Lya, with nearly all of that arising from Source A. 
Second, in Section 3.2 we argued that the metal line emis- 
sion is not cospatial with the Lya+He ii nebula, and that 
the observed metal lines likely derive from the region of 
Source A. Below we discuss the likelihood that an AGN 
at the position of Source A is powering the Lya+He ii 
nebula. 

If we assume that the Night 2 metal line emission stems 
primarily from Source A, we can compare the measured 
line ratios to those of well-studied galaxy populations. 
At face value, the C iii] and C iv emission lines associ- 
ated with the region around Source A are suggestive of an 
AGN, however we find that the line ratios are more con- 
sistent with that seen in non- AGN LBGs (jShaplev et all 
12003( 1. The measured ratios of C iv/LyQ;=0.05±0.03 and 
C ili]/Lya=0.11±0.02 (uncorrected for Lya absorption) 
are likely upper limits due to the uncertainties in the 
Lya measurement. Even so, the C iv/Lya ratio is a 
factor of 4-5 lower than what is seen in LBGs flagged 
as n arrow-line AGN (C iv/LyawO.25; IShaplev et all 
2003f) and in local Seyfert galaxies (C iv/Lya«0.21; 



Ferland fc Osterbrockl [19861 ). In addition, if Lya is un- 



derestimated by a factor of two, the corrected C iv/Lya 
and C III] /Lya ratios would match those of non- AGN 
LBGs. Furthermore, the ratio of C iii]/C iv=:2.3±1.2 
is in agreement with that found for non- AGN LBGs 
(C iii]/C iv«2.5; IShaplev et all l2003f ) and is a factor 
of four higher than expected from an narrow-line AGN 
spectrum (C iii]/C ivw0.05: IShaplev et al.ll2003f ). sug- 
gesting Source A has a softer ionizing continuum. 

There is currently no evidence from the multi- 
wavelength SED that Source A is an AGN. Existing X- 
ray coverage of the field reveals no X-ray detection at 
the position of the system {Fx > 1.5 x 10^^^ erg s^^ 
cm ~^ or Lx > 2 .7 x 10'''' erg s"\ 2-7 keV observed; 



iKenter et al.l 120051 ) but is too shallow to rule out lower 
lumiosity Seyfert galaxies. The typical luminosities of 
Seyfert galaxies: L^ya = 10^^ — 2 x 10'*'' erg s~^ and Lx 
(0.5-4.5 keV) = 5 x lO-^^ - lO^^ erg/s, with LLya/Lx 
ratios of - 0.1-2 (1 Crisi[l98l . If we combine the 
Lbya/Lx ratio as an upper limit (since the ratio for an 
extended nebula will be smaller due to the smaller cov- 
ering fraction of the gas) with our measured Lya lumi- 
nosity (5.4 X 10''^ erg s"-'), we should expect Lx (0.5-4.5 
keV)> 2.7 X 10'*^ - 5.4 X 10'*'' erg s"\ which is at or below 
our current X-ray detection threshold. 

The optical and MIR photometric measurements for 
Source A are listed in Tabled We measured the broad- 
band optical fluxes of Source A from NDWFS using 2.1" 
diameter apertures and applied aperture corrections of 
[1.06, 1.40, 1.07] in the Bw, R, and / bands. The Source 
A IRAC photometry comes from the Spitzer Deep Wide- 
Field Survey (SDWFS; 3.5" diameter aperture with point 
source aperture corrections of [1.38,1.38,1.38,1.42]). The 
IRAC colors of Sour ce A lie outs i de the AGN color-color 
select ion regions of iStern et all (|2005D and iLacv et all 
12004) , and the probability t hat an X-ray AGN w ill have 
these IRAC colors is smafl (|Goriian et all I2008D . The 



IRAC colors of Source A are [3.6] - [4.5] = 0.52 ± 0.03 
(Vega) and [5.8] - [8.01 = 0.131 ±0.11 ( Vega) (in "Region 
B" of Figure 4 of lCoriian et al.ll2008t ). The percentage 
of X-ray sources down to the XBootes limit with IRAC 
colors in this region is 4%. In a small portion of the 
survey with 10 times deeper X-ray coverage, the distri- 
bution of sources in IRAC color-color space is similar, 
i.e., there are ver y few X-ray sources with the IRAC col- 
ors of Source A (|Goriian et al.|[2008l ). In contrast, these 
IRAC colors ar e consistent w ith star-forming galaxies at 
z = 1.25-1.75 (jPonlev et a l. L200 8). the redshift range of 
our source. Similarly, the IRAC SED of Source A shows 
the 1.6/im bu mp rather than the power-la w typical of ob- 
scur ed AGN (|Ai onso-Herrero et al]|2006[ ). indicating the 
MIR SED is dominated by stellar emission. All of the 
other sources in the vicinity show similar non-power-law 
SEDs. 

PRGl is also undetected at longer wavelengths: the 
MIPS 24 /im non-detection corresponds to an upper limit 
of 51 /iJy (Icr), and the Westerbork 20 cm survey non- 
detection yields a 5 a limit at 3 GHz in the restframe of 
6.4 X 10^3 W Hz-i (Ide Vries et al.l[200l . weU below the 
realm of high redsh ift radio galaxies (~ 10^^ W Hz~*; 
ISevmour et"ani2007f) . 

On the other hand, the energetics of the nebula sug- 
gest that an AGN at the projected position of Source 
A could explain the observed continuum emission if the 
AGN is highly obscured to our line-of-sight. If we make 
the assumption that an AGN at Source A is power- 
ing the He ii emission, we can estimate the amount 
of Bw continuum emission we expect from the neb- 
ula due to illuminatio n by the AGN. Scaling a stan- 
dard AGN template (jElvis et al.l I1994D to match the 
ionizing photon flux implied by the He ii (Q(He+) w 
6.2 X 10^3 photons s~*), we estimate that the Bw lu- 
minosity from the AGN striking the cloud should be 
LBwMcident ^1.9x10"* erg s"* A~*. This incident AGN 
continuum flux — the maximum possible contribution 
from the AGN to the observed continuum emission — 
is roughly the same as the measured Bw continuum 
emission from the nebula {LB„,nebuia sal.lxlO'** erg s^* 
A"*). In reality, the continuum observed within the neb- 
ula will arise solely from two-photon (27) continuum and 
scattering of AGN light with no contribution directly 
from the AGN itself. Given the predicted He 11 equiv- 
alent width from our Cloudy models of AGN illumina- 
tion, we estimate the expected Bw two-photon contin- 
uum within the nebula to be Lb„,2-i ~2.0x10''° erg s~*, 
roughly 18% of the observed Bw continuum of the neb- 
ula. Electron scattering of AGN light is expected to have 
a scattering optical depth of t = Uelar— 0.1-1.1, where 
TT-e wl-lO cm~'^ and I ~ 56 kpc is the approximate path 
through the nebula, corresponding to a Thompson scat- 
tering continuum of Lb^v , Thompson <1.9xl0''* crg s"* 
A"* X (1 - e~^) = 1.3x10"* erg s~* A"*, comparable 
to the observed Bw continuum emission. Although dust 
scattering, which is much more efficient than Thomp- 
son scattering, may contribute as well, large quantities 
of dust would be difficult to reconcile with the large ob- 
served Lya and He 11 luminosities. Questions remain, 
but given the limitations of our discovery data it ap- 
pears that an AGN at the position of Source A that 
is unobscured along the line-of-sight to the nebula, but 
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completely hidden from our viewpoint, could explain the 
observed Bw continuum emission. Correcting for the 
maximal (since we are assuming the projected separa- 
tion) solid angle subtended by the cloud as seen from the 
Source A {dQ « 3.1 steradians), wc find that the mini- 
mum intrinsic AGN luminosity needed to produce this in- 
cident flux is LB„,min{AGN) >1.9xl0'*i crg s^'^xiTT/dn 
erg s~^ w7.6xl0^^ erg s~^x A~^. Source A has 

a B\Y luminosity of 1.8x10^° erg s~^x A~^, which is 
over an order of magnitude fainter than this firm lower 
limit. The AGN would therefore need to be very highly 
obscured along our line-of-sight to match the existing ob- 
servational constraints. 

In conclusion, an AGN at the position of Source A is a 
plausible source of the hard ionizing radiation. However, 
there is no visible AGN in the vicinity of the nebula, so 
it would have to be highly obscured along our line-of- 
sight. An AGN in the vicinity of the nebula must be in 
a radio-quiet phase and so highly obscured that even the 
observed MIR SED is dominated by light from the host 
galaxy. Deep optical and near-infrared spatially-resolved 
spectroscopy will be required to resolve this issue. 

4.4.2. Shocks 

The observed ratios are inconsistent with shock ion- 
ization in solar metallicity gas, which typically produces 
much higher C ill] /He ii and C iv/C III] ratios (e.g., 
^3-25 and ^1-10 for shock velocities of 500-150 km s~^; 
iDopita fc Sutherlandlll996l ) along with strong N v]A1240. 
Stronger shock velocities are inconsistent with the nar- 
row width of the He ii line {vfwhm ^500 km s~^). 

4.4.3. Wolf-Rayet Stars 

Strong He ii emission is seen in both of the major 
classes of Wolf-Rayet stars (WN and WC), but it is ac- 
companied by strong Si iv, N v], and/or other metal 
emission lines. WN stars show strong Si ivAA1393,1402 
relative to He ii, C iv, and C ill], which we can rule 
out with our discovery spectra. WC stars usually show 
C iv/C III] ratios greater than 1, medium-strength Si iv, 
and a large number of other metal lines (e.g., Fe), all of 
which are inconsistent with our observations. The spec- 
trum of a typical 'W-R galaxy' (a galaxy with spectrum 
dominated by Wolf-Rayet features) effectively averages 
over hundreds or thousands of W-R stars, but none-the- 
less typically shows Si iv, C iv, He li, and N v] emis- 
sion with P-Cygni profiles due to the e ffects of strong 
stellar winds (e.g., iLeitherer et al.lfr996l ). The compos- 
ite Lyman Break Galaxy (LBG) spectrum measured at 
z ^ 3 also shows Si IV and C iv with P-Cygni profiles 
from the stellar winds of massive stars, as well as broad 
He II {vfwhm ~1500 km s~^), which the authors argue 
is mo st likely due to the f ast, dense winds of Wolf-Rayet 
stars (iShaplev et al.l[2003l) . In our case, the narrow width 
of the He ii line {vfwhm ^500 km s~^) and the absence 
of other important Wolf-Rayet features (e.g.. Si iv, N v]) 
rules out the idea that the He ii emission is coming from 
Wolf-Rayet stars. 

4.4.4. Population II Star Formation 

The existence of diffuse, blue continuum that is cospa- 
tial with the Lya and He ii line emission (Figures [Ij 
is suggestive of a distributed source of ionization such 



as spatially-extended star formation. However, the ob- 
served He il/Lya line ratio is inconsistent with ioniza- 
tion by a standard Population II SED (e.g., Z = 0.001 — 
I/2OZ0). The expected Q(He+)/Q(H) ratio for a Pop II 
model is several orders of magnitude lower than observed 
(Schaerer. .2003.) . Predictions from our Cloudy models 
indicate that the dearth of hard ionizing photons trans- 
lates into negligible He 11 emission and He ii/Lya ratios 
that are a factor of 1000 lower than observed. Further- 
more, the large observed He 11 equivalent width can only 
be produced by very low metallicity stellar populations 
(lSchaerei]l2n03l i. 

4.4.5. Population III Star Formation 

The only way to explain strong He 11 emission with 
ionization from a stellar population is to invoke very 
young ages and very low metallicities. The rest-frame 
equivalent width of Lya (EWw294A) is higher than or 
comparable to what is expected from solar me tallicity 
and Pop II stellar populations (jSchaereil l2003f l. but it 
is by no means the largest obse rved in Lyg-em itting 
galaxy surveys (EWmax >300A: [Pawson et ahl [2007). 
The measured equivalent width of He 11, however, is 
very large (EW«37A) in the context of stellar popula- 
tions and only consistent with the youngest (<2 Myr) 
and lowest metallicity stars (<10 ~^ Zq), in the ab- 
sence of mass loss dScha ereil l2003l ) . Here we chose a 
Pop III model with a moderately top-heavy IMF from 
iTumlinson et al.l (|2006D that has a peak at 10 Mq (their 
case A); a more top-heavy IMF will tend to boost the 
He ii/Lyck ratio due to the additional hard ionizing pho- 
tons. The ionizing photon fiux implied by the Lya and 
He II measurements from PRGl implies a Pop III cluster 
mass of 13-100 xlO^ M0 and a Bw continuum flux of 
LBw,tncident «1.7xl0''i erg s^^ A"\ roughly equal to 
the observed continuum {LBw,nebuia «l.lxlO''^ erg s~^ 
A~^). In comparing with our Cloudy models, how- 
ever, we find that the observed He ii/Lya ratio (RiO.12) 
is much higher than expected from our Pop III model 
(«0.004). However, it is important to keep in mind that 
model uncertainties, such as the assumed IMF and the 
effects of mass loss, could have a large effect on the pre- 
dicted line luminosities. In addition, the possible under- 
estimation of the Lya flux discussed in Section 4.1 as 
well as a density-bounded geometry will tend to boost 
the observed He ii/Lya ratio. 

On the face of it, it would be surprising to find such low 
metallicity star formation at such a low redshift. How- 
ever, we cannot conclusively rule out the Pop III scenario 
on this basis alone. While the mean metallicity of the 
Universe increases with time, several theoretical models 
of Pop III star formation have suggested that significant 
metall icity inhomogeneities may exist even at low red- 
shifts (iTornatore et al.ll2007f l. These models predict that 
while the metallicity is quickly enriched at the center of 
collapsed structures, low metallicity regions can persist 
on the periphery over longer timescales, allowing Pop III 
star fori nation to proceed wel l after the epoch of reion - 
ization (jTornatore et al.ll2007l: IScannapieco et al.ll2003D . 
At z « 2.3 (ro ughly the redshift mid point of our Lya 
nebula survey) , ITornatore et al.l (|2007 ') predict a Pop III 
SFR density of 1.3x10"'' Mq yr"^ Mpc"^. If we were to 
assume that this Lya-|-He 11 nebula is powered by very 
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low metallicity star formation, use the Lya line emission 
(extrapolated to the entire nebula and corrected for ab- 
sorption) to estimate a SFR (SFR = LLya / 1.26x10^^ 
5.4x 10-^3 / 1.26x10^^ erg s-i « 42.9 MQyr-i;[Kmnicuti 
Il998f l. and ignore any incompleteness of our survey, we 
would calculate that one nebula of this kind within our 
survey volume (1.3x10* Mpc"^), corresponds to a Pop III 
SFR density of ~ 3.3 x 10"'^ Mq yr'^ Mpc'^. This very 
rough estimate based on a single source is of the same 
order of magnitude as the model predictions (within a 
factor of 3) despite the large theoretical uncertainties in 
the conditions regulating Pop III star formation and IGM 
enrichment at all redshifts. 

In reality, the amount of low metallicity gas and Pop III 
star formation as a function of rcdshift will depend on 
the efficiency of star formation in different environments 
and the efficiency with which stars pollute their environ- 
ment, i.e., the feedback efficiency. A more realistic treat- 
ment of galactic winds in cosmological simulations (e .g., 
momentum-driven winds; Oppcnhcimcr & Dave 2006!) is 
shown to suppress the metallicity in all gas phases rel- 
ati ve to the consta n t win d scenario similar to that used 
in iTornatore et al.l ()2007l ). which may in turn increase 
the fraction of Pop III star formation at any given red- 
shift. On the other hand, only a few stellar genera- 
tions would be required to produce the mass of carbon 
in the nebula. Assuming the derived ionized gas mass 
(~ 8 X 10* M0), the simulated yields for Pop III super- 
novae (^0.2-1 Mp, C p er 20-40 Mq supernova progenitor; 
iTominaga et al.l[2007f l and instantaneous mixing, enrich- 
ing the cloud to z ^ 10~^Zq would only require of order 
750 supernovae (2 — 4 x lO'' Mq). This is several or- 
ders of magnitude lower than the fraction of the cluster 
mass implied by the Lya and He II luminosities that is in 
stars with M > 20Mq for the assumed top-heavy IMF: 
/(M > 20Mq) = 0.26 or 3-27x IO^Mq. Thus it is likely 
the nebula would be polluted to the observed levels over 
a relatively short timescale. 

Large uncertainties remain in our understanding of 
Pop III star formation. The expected fraction of Pop III 
star formation at a give epoch is unclear, feedback effi- 
ciency of these first stars is largely unconstrained, and 
there are a wide range of possible Lya/He 11 ratios, de- 
pending on stellar mass loss, star formation history, and 
the complicated effects of Lya radiative transfer. Despite 
the strong Lya-l-IIe 11 signature in this source, we cannot 
conclusively confirm or rule out Pop III star formation 
as a source of ionization with the current observations. 

4.4.6. Gravitational Cooling Radiation 

A number of theoretical papers have addressed the 
issue of gravitational cooling radiation, i.e., the cool- 
ing of gas as it collapses within the dark matter po- 
tential and heats to T « lO'' K. Metal-line cooling 
is possible for gas with T < 10^ K, but it is ~1000 
times less efficient than cooling via L ya emission, un- 
less Z > 0.1 Zq (jHaiman et all |2000D . Thus for low 
metallicity gas, str ong Lya emission is expected to dom- 
inate the cooling (jFardal et al.l [2OOII ). He 11 emission 
may be important as well, but its contribution relative 
to Lya is unclear due to uncertainties in how much Lya 
will escape the system. From a semi-analytic calcula - 
tion assuming monolithic collapse, iHaiman et al.l (120001 ) 
suggested that a cloud that is sufficiently metal-poor 



will radiate 10% of the energy via He iiA304A; this 
corresponds t o a He llA1 6 40/Ly a of ^0.01, where we 
have followed lYang et al] (|2006[ ) and adopted a ratio 
of He iiA1640/He iiA304p a 0.10, i.e., the case B val- 
ues of IStorev fc Hummeil (jl995) extrapolated to the 
low density limit. Using Parallel TreeSPH simulations, 
lYang et al.l (|2006D predicted ratios of He iiA1640/Lya< 
0.10 — 0.001 for gravitationally coofing clouds, depending 
on the degree to which self-shielding of the gas reduces 
the Lya flux. This range approaches the observed ra- 
tio for PRGl, but due to subsequent corrections to the 
treatment of star-for ming gas in more recent simulations 
the predictions from lYang et al.l (120061 ) are likely overes- 
timates (Y. Yang 2008, private communication). Thus, 
the observed He ii/Lya ratio appears to be higher than 
predicted for gravitational cooling radiation, but again 
the uncertainties in both the theoretical predictions and 
in our Lya measurement make it difficult to draw firm 
conclusions. 

Some of the theoretical predictions for gravitation- 
ally cooling clouds are consistent with our observations. 
The observed He 11 line width {a^ < 250 km s^^) 
is more consistent wit h gravitational co oling radiation 
{ay < 400 km s'^ e.g.. lYang et al.ll2006[ ) than with out- 
flows. In addition, theoretical simulations of the red- 
shift evolution predict a peak in the number density of 
gravitational cooling Lya nebulae at z « 2, consistent 
with o ur discovery of a Lya-I-He 11 nebula at z « 1.7. 
While iDiikstral (|2007D predicts rest-frame Lya equiva- 
lent widths of >1000A for gravitationally cooling clouds, 
he notes that the observed rest-frame equivalent width is 
likely to be reduced by a factor of 5-10 due to scattering 
of Lya photons in the intergalactic medium, more in line 
with our measurements. 

However, a number of other theoretical predictions of 
gravitational cooling radiation do not fit our observa- 
tions. First, the simulations do not reproduc e the relative 
sizes o f the observed Lya and He 11 regions. iFardal et al.l 
(|2001l ) could not reproduce the sizes of the largest ob- 
served Lya nebulae unless they invoked resonant scat- 
tering of the Lya emission. In that case, the He 11 emis- 
sion shoul d be more centrall y-concentrated than the Lya. 
Similarly, Yang et all ()2OO60 suggested that He 11 regions 
would likely be too small to resolve with current observa- 
tional facilities. In contrast, in our Lya-I-He 11 nebula the 
He II region appears to be comparable in size to the Lya 
region (w45 kpc). The theoretical papers also predict 
that Lya nebulae will only be present as a halo around 
a ma ssive galaxy (|Fardal et al.l 120011 : iFurlanetto et al.l 
[2001 . While there are a few faint sources around the 
edge of our nebula, and there may be unresolved low 
surface brightness clumps within the cloud, the multi- 
wavelength imaging shows no evidence for a central mas- 
sive galaxy in this system. 

Given the uncertainties in the theoretical predictions 
and the limits of our discovery data it is difficult to as- 
sess the applicability of the gravitational cooling model. 
The weight of the current evidence disfavors gravitational 
cooling as the sole explanation for the line emission, but 
it is certainly possible that the nebula is powered by 
multiple processes, with gravitational cooling radiation 
playing a larger role on the outskirts and photoioniza- 
tion from stars or AGN dominating the ionization of the 
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central regions. 

5. IMPLICATIONS 

Strong Lya+He II in the absence of strong metal lines 
has been championed as a unique observational signa- 
ture of primordial galaxy formation (e.g., Pop III star 
formation or gravitationally cooling clouds) , but the dis- 
covery of this Lya-I-He li nebula suggests that the situ- 
ation can be much more complex in actual astrophysi- 
cal sources. Occam's razor suggests that the most likely 
power source is an AGN at the position of Source A. An 
analysis of the existing data shows no obvious evidence 
of a powerful AGN in the vicinity, so to explain the ion- 
ization of the nebula, an AGN would need to be highly 
obscured even in the mid-infrared. The line ratios rule 
out ionization by shocks, Wolf-Rayet stars, and Pop II 
star formation. Pop III star formation would provide 
the necessary hard ionizing continuum to explain the ob- 
served line ratios along with a natural explanation for 
the spatially-extended continuum emission, but this sce- 
nario appears unlikely given the low redshift. Despite the 
compelling Lya-I-He ii signature, the complexity of this 
source and the large uncertainties in theoretical predic- 
tions preclude a more definitive conclusion. The contri- 
bution from gravitational cooling radiation is similarly 
unclear, although the morphology of the nebula (with 
no central compact source) and relative sizes of the Lya 
and He ii emitting regions disfavor this hypothesis as a 
dominant mechanism. 

One of the most important implications of this discov- 
ery is that it demonstrates we must be careful when using 
strong Lya-I-He ii emission as a tracer of Pop III star for- 
mation. Surveys looking specifically for this Lya-I-He ii 
signature are ongoing. While the low redshift of PRGl 
allows for the detailed multi-wavelength follow-up nec- 
essary in order to understand the power source and the 
metallicity of the gas, the same cannot be said for on- 
going Lya-|-He ii emitter surveys at higher redshift that 
will lack longer wavelength coverage and be sensitivity- 
challenged. It is extremely important to note that if this 
Lya-I-He ii had been discovered at high redshift, it would 
have been easy to jump prematurely to the Pop HI con- 
clusion. The discovery of a Lya-I-He ii nebula at z w 1.67 
therefore provides an ideal opportunity to evaluate the 
extent to which strong Lya-|-He ii emission can be used 
as a unique tracer of Pop HI star formation and under- 
scores the importance of using care when interpreting a 
strong Lya-|-He II signature in the absence of more ex- 
tensive multi-wavelength data. 

6. SUMMARY 

We have discovered a Lya nebula at z « 1.67 (the low- 
est redshift Lya nebula known) that shows strong, diffuse 
He II emission and weak/negligible C ill] and C iv emis- 
sion. From the line ratios we derive evidence that this 



nebula may contain low metallicity {Z < 10^^ — lO^'^Z©) 
gas, depending on the unknown ionization parameter, 
that is being illuminated by a hard ionizing continuum, 
either due an AGN or a very low metallicity stellar pop- 
ulation (Pop III), by gravitational cooling, or some com- 
bination thereof. The softer continua of Population I 
and II stars can be conclusively ruled out along with 
ionization by shocks and Wolf-Rayet stars. No obvious, 
unobscured, powerful AGN is seen in the system; thus if 
an AGN is responsible, it must be highly obscured along 
our line-of-sight. Despite the strong Lya-I-He ii signa- 
ture, our detailed analysis of the discovery data shows 
that dedicated follow-up observations will be required 
in order to draw firm conclusions about the dominant 
source of ionization for the nebula and better constrain 
the metallicity. This is the first time that the predicted 
observational signature of very low metallicity (Pop III) 
star formation — strong Lya and He ii in the absence of 
strong metal lines — has been seen in a spatially-extended 
source; however, the complex nature of the nebula and 
the fact that such complexity becomes increasingly diffi- 
cult to discern with redshift suggest that studies at high 
redshift will need to use caution when interpreting future 
Lya-I-He il discoveries. 



We are grateful to J. Tumlinson for providing the 
Pop III evolving spectra, to E. Le Floc'h and the IRS 
GTO and MIPS GTO teams for providing the MIPS 
data, to Gary Ferland for useful advice on the photoion- 
ization modeling, and to K. Finlator for observing as- 
sistance and many helpful discussions. D. Stern, M. 
Ashby, M. Brodwin, and the rest of the SDWFS team 
are thanked for access to the most current Spitzer IRAC 
imaging of the Bootes field. We also thank the anony- 
mous referee for suggestions that improved the clarity of 
this paper. We would like to acknowledge the expert as- 
sistance of the staff of the MMT Observatory (especially 
John McAfee, Ale Milone, and G. Grant WiUiams). M. 
P. was supported by an NSF Graduate Research Fellow- 
ship and a P.E.O Scholar Award. This research builds on 
data from the NOAO Deep Wide-Field Survey (NDWFS) 
as distributed by the NOAO Science Archive. NOAO is 
operated by the Association of Universities for Research 
in Astronomy (AURA), Inc. under a cooperative agree- 
ment with the National Science Foundation. Some of 
the data presented in this paper were obtained from the 
Multimission Archive at the Space Telescope Science In- 
stitute (MAST). STScI is operated by the Association 
of Universities for Research in Astronomy, Inc., under 
NASA contract NAS5-26555. Support for MAST for 
non-HST data is provided by the NASA Office of Space 
Science via grant NAG5-7584 and by other grants and 
contracts. 



REFERENCES 



Alonso-Herrero, A., et al. 2006, ApJ, 640, 167 

Bromm, V., Kudritzki, R. P., & Loeb, A. 2001, ApJ, 552, 464 

Dawson, S., et al. 2004, ApJ, 617, 707 

Dawson, S., Rhoads, J. E., Malhotra, S., Stern, D., Wang, J., Dey, 
A., Spinrad, H., & Jannuzi, B. T. 2007, ApJ, 671, 1227 

Dey, A., van Breugel, W., Vacca, W. D., & Antonucci, R. 1997, 
ApJ, 490, 698 



Dey A., et al. 2005, ApJ, 629, 654 

Dijkstra, M. 2007, ArXiv e-prints, 711. l arXiv:0711.2698l 

di Serego Alighieri, S., Kurk, J., Ciardi, B., Cimatti, A., Daddi, E., 

& Ferrara, A. 2008, ArXiv e-prints, 807, arXiv:0807.4634 
Donley J. L., Rieke, G. H., Perez-Gonzalez, P. G., & Barro, G. 

2008, ArXiv e-prints, 806, arXiv:0806.4610 
Dopita, M. A., & Sutherland, R. S. 1996, ApJS, 102, 161 



A Lya+Heii Nebula at z w 1.67 11 



TABLE 1 

PRGl Spectroscopic Measurements 



LyaA1216 NvA1240 SiivA1398 CivA1549 HeiiA1640 Cin]A1909 NcivA2424 
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Redshift 
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FWHM<,b, [A] 
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TABLE 2 
PRGl Photometric Measurements 





iNeDuia ' 
10-30 erg s-l cm-2 Hz"! 


Source A'^ 
10-30 erg s-l cm-2 Hz-l 




3.58 ± 0.24 


2.07 ± 0.14 


R 


4.03 ± 0.49 


4.60 ± 0.33 


I 


6.62 ± 0.88 


8.35 ± 0.47 


3.6/im 




466.56 ± 34.13 


4.5/im 




488.70 ± 32.40 


5.8/im 




329.67 ± 31.40 






194.13 ± 25.60 



^ Flux measured within the same 1.5 X 5.0" aperture as the spectroscopic 
measurements.^ The contribution of line emission for Night 1 (Night 2): 
He II contributes 8% (7%) and C IV < 2% (3%) of the Bw emission, 
and Ne iv] contributes < 2% (6.5%) of the J?-band emission. Wc sec no 
contribution from line emission out to the middle of the /-band (8300A, 
the extent of our spectroscopic coverage).'^ Optical Bw , R- and / fluxes 
were measured within 2.1 ' diameter apertures. Mid-infrared fluxes (3.6, 
4.5, 5.8, 8.0 fj,m) were measured within 3.5 diameter apertures. Aper- 
ture corrections are discussed in the text. 



TABLE 3 
Emission Line Ratio Comparison 





Lya/Heli'' 


Civ/Hell 


Cm] /Hen 


Civ/Ciil] 


Reference 


Radio Galaxy Halos (Composite) 


11.7 


1.75 


0.73 


2.4 


Humphrey et al. 2008*= 


Radio Galaxy Halos (Mean) 


9.80 ± 5.69 


1.59 ± 0.56 


1.06 ± 1.05 


2.22 ± 1.17 


Humphrey ct al. 2008 


Stacked Lya Blobs at 2 3 


11.11 ± 9.88 








Saito et al. 2008 


Lya Blob at 2 2.7 


7.62 ± 0.08 


1.02 ± 0.01 


0.12 ± 0.02 


8.34 ± 1.67 


Dey et al. 2005 


PRGl Night 1 


8.00 ± 2.32 


< 0.22 


< 0.19 




This study 


PRGl Night 2 


7.59 ± 1.43 


0.36 ± 0.20 


0.83 ± 0.20 


0.44 ± 0.24 


This study 



^ No correction has been applied for Lya absorption.^ Errors on line ratios from composite spectrum not given. la upper 
limits. 
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Fig. 1.— GALEX (FUV and NUV), NDWFS broad-band optical (Bw, R, and /), and IRAC {3.6ftm, 4.5;tm, 5.8;jm, and 8.0/im) 
postage stamps for PRGl. Images are all 1' on a side and centered on the coordinate location 14:35:12.4.39 +3-5:11:07.16 (J2000). The 
spectroscopic slit is shown with a position angle of 81.2°. PRGl was selected as a Lya nebula candidate due to the diffuse blue emission 
in the imaging. However, in this case Lyo is in fact outside the B^ band; the B^^ flux is instead dominated by diffuse continuum 
and He ii emission. The diffuse blue continuum emission is clearly visible in the NDWFS Biyband imaging. Several compact sources in 
the near vicinity of the nebula have IRAC counterparts, labeled A, B, and C. 
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Fig. 2. — Night 1 (UT 2008 June 08) 2D and ID spectra showing strong Lya and He II emission from a single half hour exposure 
(1.5x5.0" aperture). The spectra have been smoothed by 0.84" spatially and by 5.8 A in the spectra dimension. The filter curve is the 
By^ bandpass convolved with the atmospheric transmission (blue). The la error spectrum is overlaid (red). 
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Fig. 3. — Night 2 (UT 2008 June 09) 2D and ID spectra showing strong Lyo and He II emission, weak C III], and marginally-detectable 
C IV from a combined 1.5 hour exposure (1.5x5.0" aperture). The spectra are smoothed by 0.84" spatially and by 5.8 A in the spectral 
dimension. The filter curve is the Bw bandpass convolved with the atmospheric transmission (blue). The la error spectrum is overlaid 
(red). When comparing to Figure [2] note that a 180° change in the slit position angle between the two nights caused a flip in the angle of 
the spectral trace across the detector. 
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Night 1 Night 2 
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Position (") Position (") 

Fig. 4. — Lya and He II spatial profiles with errors as a function of position along the slit from the Night 1 (UT 2008 June 08) and Night 
2 (UT 2008 June 09). As C in] was not detected on Night 1, the C III] profile is shown for Night 2 only. He II and C III] profiles scaled by 
a factor of 6 are also overplotted to allow easier comparison with the Lya profile. The He II region appears to be comparable in size to the 
Lycf region. The shaded area represents the approximate position of Source A, a source that was off-slit but that may have contaminated 
the Night 2 observations due to poor seeing and windy conditions. The discrepancy between the spatial profiles and the offset of the C III] 
spatial profile relative to the hya are both consistent with the hypothesis that the Night 2 spectrum is contaminated by light from Source 
A and that Source A may be the primary source of the metal line emission. 
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Fig. 5. — Lyo (black solid line) and He II (blue dashed line) line profiles as a function of observed wavelength centered on the systemic 
redshift of the system, as measured from the He II line. The observed Lya is shown with the mirror image of the long-wavelength half 
of the line superimposed (black dotted line). A Gaussian fit (red solid line) indicates that the Lya line may be absorbed by ~41%. The 
corresponding velocity offsets for Lya are given on the top axis. 
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Fig. 6. — Velocity and velocity dispersion spatial profiles of the Lya line for Night 1 (UT 2008 June 08; black filled circles) and Night 2 
(UT 2008 June 09; red open circles) derived from spectra extractions taken in 2 pixel (0.56") spatial bins and corrected for the instrumental 
resolution. The velocity zeropoint was set using the redshift calculated from the He II line centroid in the full 5.0" aperture extraction. 
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Fig. 7. — C Ill]/He II versus C iv/He II line ratios for a grid of Cloudy models over a range of SEDs, ionization parameters, and gas 
metallicities. The models represent a cloud of gas with LogZ/Zp) =[-3, -2, -1, -0.5, 0] being ionized with an ionization parameter of 
LogU =[-3, -1, 0] by one of three SEDs: an AGN SED (left), a Tumlinson et al. (2006) Population III galaxy model SED (middle; Z=0, 
1 Myr burst, moder ately top-heavy IMF d iscussed in the text), and a Starburst99 Population II galaxy model SED (right; Z=0.001, Salpeter 
IMF, 1 Myr burst; ILeitherer et al.lll996l) . The Icr limits on the line ratios of the nebulosity are shown for Night 1 (solid line) and Night 2 
(dashed line). 



